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ABSTRACT 
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■ Aims. Observations of Kepler's supernova remnant (G4.5+6.8) with the H.E.S.S. telescope array in 2004 and 2005 with a total live time of 13 h 
' are presented. 

Methods. Stereoscopic imaging of Cherenkov radiation from extensive air showers is used to reconstruct the energy and direction of the incident 
OO gamma rays. 

• Results. No evidence for a very high energy (VHE: >100 GeV) gamma-ray signal from the direction of the remnant is found. An upper limit (99% 
!L/ " confidence level) on the energy flux in the range 230 GeV - 12.8 TeV of 8.6 x 10~ 13 erg cirT 2 s is obtained. 

5^ ' Conclusions. In the context of an existing theoretical model for the remnant, the lack of a detectable gamma-ray flux implies a distance of at least 
, 6.4 kpc. A corresponding upper limit for the density of the ambient matter of 0.7 crrr 3 is derived. With this distance limit, and assuming a spectral 
j_j ■ index T = 2, the total energy in accelerated protons is limited to E p < 8.6 x 10 49 erg. In the synchrotron/inverse Compton framework, extrapolating 

■ the power law measured by RXTE between 10 and 20keV down in energy, the predicted gamma-ray flux from inverse Compton scattering is 
below the measured upper limit for magnetic field values greater than 52 pG. 

Key words, gamma rays: observations - ISM: supernova remnants - ISM: individual objects: (Kepler's SNR, SN1604, G4.5+6.8) 



1. Introduction 

It is widely believed that the bulk of the Galactic cosmic rays 
(CR) with energies up to at least several 100 TeV originates from 
supernova explosions (see for example Drury et al. 1994). This 
implies copious amounts of very high energy (VHE: >100 GeV) 
nuclei and electrons in the shells of supernova remnants (SNRs). 
These particles can produce VHE gamma rays in interactions of 
nucleonic cosmic rays with ambient matter, via inverse Compton 
(IC) scattering of VHE electrons off ambient photons, as well 
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as from electron Bremsstrahlung on ambient matter. Therefore 
SNRs are promising targets for observations of VHE gamma 
rays. 

In October 1604 several astronomers, among them Johannes 
Kepler, observed a "new star" which today is believed to have 
been a bright supernova (SN) at the Galactic coordinates / = 
4.5° and b = 6.8°. The remnant of this supernova has since 
been a target of observations covering the entire electromagnetic 
spectrum. In the radio regime, [Dickel et al. (1988)] determined a 
mean angular size of ~ 200" and a mean expansion law R oc f 50 , 
where R is the radius and t is the time. However, the expan- 
sion parameter x = Rt/R varies considerably around the SNR 
shell, 0.35 < x < 0.65, possibly indicating spatial inhomogen- 
ities in the circumstellar gas density. In a very recent paper by 
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|Vink (2008)| these properties, and the general asymmetry of the 
remnant, have been basically confirmed through X-ray measure- 
ments. They also allowed the analysis of a high-velocity syn- 
chrotron filament in the eastern part of the remnant with x = 0.7. 
In addition, the distance d to the SNR is still under de- 



bate. Reynoso & Goss (1999) report on an HI absorption fea- 
ture in VLA data and use the Galactic rotation model of 
|Fich et al. (1989)| to calculate a lower limit d > (4.8 ± 1.4)kpc. 
They also give an upper limit on the distance due to the lack 
of absorption by an HI cloud at 6.4 kpc. The authors remark 
that these values involve uncertainties because of the prox- 
imity of Kepler's SNR to the Galactic center. In contrast, 
|Sankrit et al. (2005)| and subsequently |Blair et al. (2007)| have 
given a lower source distance of d — 3.9(+1.9 - 0.9) kpc, from 
an absolute shock velocity ~ 1660 + 120kms _1 derived from 
the Ha emission line width of a Balmer-dominated filament that 
is located in the northwestern region. The line broadening, taken 
as an indication of the downstream thermal gas temperature, was 
used to determine the shock velocity. We shall return to this 
question in the discussion section. 

Finally, the type of the supernova is not undisputed. From 
the reconstructe d light curve [ Baade (1943) claimed that it was a 
type la SN, but |Doggett & Branch (1985)| argued that the li ght 
curve is also consistent with a type II-L. |Smith et al. (198"9)| and 



Kinugasa & Tsunemi (1999)|observed a relative overabundance 



of heavy elements that agrees with type la nucleosynthesis mod- 
els, while |Decourchelle & Ballet (1994)| saw more evidence that 
Kepler's SNR is the remnant of a core-collapse SN. Its posi- 
tion, 500 - 750 pc above the Galactic plane, is more consis- 
tent with a type la than a type II SN, as a SN of the latter 
type is expected to be confined to the region of high gas den- 
sity found in the plane. However, in the case of a core-collapse 
event this might be explained through the model of a runaway 
star, as proposed by Bandier a (1987)[ More recently, theoreti- 
cal modeling of the detailed thermal line spectra obtained with 
XMM dCassam-Chenai et al. 2004i i led |Badenes et al. (2005)| to 
the conclusion that the X-ray sp ectrum is bes t fit by a type 
la SN, a view also expressed by |Blair (20 05) Most recently 



Reynolds et al. (2007) reported on deep Chandra observations 
and argued from the high abundance of iron and the very low 
abundance of oxygen that the progenitor of Kepler's SNR has 
been a type la SN. Therefore it appears that the observational 
evidence is finally converging on a type la event. 

In this paper observations of Kepler's SNR with the H.E.S.S. 
telescope array are described. An upper limit on the integrated 
energy flux above 230 GeV is derived. Combining this H.E.S.S. 
result with the theoretical predictions of |Berezhko et al. (2 006) 
suggests a lower limit on the distance, close to the upper limit 
given by |Reynoso & Goss (19 99) if Kepler's SN is a priori as- 
sumed to be of type la. 

2. H.E.S.S. Data and Analysis 

H.E.S.S. is an array of four imaging atmospheric Cherenkov 
telescopes situated in the Khomas Highland of Namibia 
dHinton 20 04). Kepler's SNR was observed with the entire tele- 
scope array between May 2004 and July 2005 for a total ob- 
servation time of 14 h. The observations were made in wob- 
ble mode, where the tracking position is offset from the source 
center (RA 17 h 30 m 42.12 s , Dec -21°28'59.9" J2000.0). Offsets 
ranging from 0.40° to 0.85° were used. The data were taken 
at zenith angles between 1° and 49°, with a mean zenith angle 
of 13°. After applying the standard H.E.S.S. data-quality crite- 
ria a total of ~13 h live time were available for the analysis. 
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Fig. 1. Left: Sky map of excess events around the position 
of Kepler's SNR with oversampling radius 0.112°; Right: 
Distribution of the squared angular distance of gamma-ray-like 
events to the center of the remnant (ON) and the center of three 
control regions (OFF) with the same distance to the pointing po- 
sition as the ON region. The vertical dotted line denotes the stan- 
dard selection cut for point sources used by H.E.S.S. 



The analysis is performed using the standard analysis techniques 
(I Aharonian et al. 20041 1 Aharonian et al. 20031 >. 

An event is counted as an ON-source event if its direc- 
tion is reconstructed within 0.112° from the direction of the 
source, given that Kepler's SNR is expected to be point-like for 
H.E.S.SQThis is a reasonable assumption as the angular size of 
the remnant in radio and X-rays wavelengths is 200"(= 0.06°). 

As the data were taken in wobble mode, the background es- 
timation can be done using OFF-source regions in the same field 
of view with the same size and offset angle (angular distance to 
the pointing position) as the source region ( |Berge et al. 2007) 1. 

A second independent analysis, used to cross-check the 
results, is based on the three-dimensional modeling of the 
Cherenkov light in the shower. (Lemoine-Goum ard et al. 2 006). 
The background estimation for this second analysis was done 
similarly. 

With the standard analysis 827 ON and 8855 OFF events 
(with a normalization of a = 0.091 1) are measured, resulting in 
an excess of 20 + 30 events. The total significance of the excess 
from the direction of Kepler's SNR (calculated using Equation 
17 of |Li&Ma(1983)| l is 0.68 standard deviations. Fig.Q]shows 
in the left panel a sky map of excess events around the posi- 
tion of Kepler's SNR and in the right panel the distribution of 
the squared angular distance of observed gamma-ray candidates 
from the center of the remnant in comparison to OFF data. The 
angular distribution of the ON events is compatible with the dis- 
tribution of the OFF events. There is no evidence for a gamma- 
ray signal from Kepler's SNR. 

The approach of Feldman & Cousins (1998) is used to cal- 
culate the upper limits on the integrated photon flux above 
230 GeV. At a confidence level of 99% an upper limit of F(> 
230 GeV) < 9.3 x 10~ 13 cirT 2 s _1 for an assumed photon in- 
dex of T = 2.0 is derived. At the same confidence level an 
upper limit on the energy flux of ^(230 GeV - 12.8 TeV) < 
8.6 x 10~ 13 ergcirT 2 s _1 in the H.E.S.S. energy range for this 
data set (230 GeV - 12.8 TeV) is derived. The assumed index of 
2 requires an upper bound for the integration range to avoid a 
divergent energy flux. 

These values depend only weakly on the assumed photon 
index for reasonable values (i.e. 2.0 < T < 3.0). 



1 The value 0.1 12° comes from a cut on the squared angular distance 
of 0.0125 deg 2 used in the standard H.E.S.S. analysis. 
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3. Discussion 

To put the observed upper limit on the gamma-ray emission 
into perspective, the above result is compared with theoretical 
expectations. Such expectations have recently been formulated 
by |Berezhko et al. (2006)| (BKV), using a non-linear kinetic the- 
ory of cosmic-ray acceleration in SNRs. This model is based 
on a time-dependent, spherically symmetric solution of the CR 
transport equation, coupled to the dynamics of the thermal gas. 
The key assumption is that the explosion was a standard type la 
event in a circumstellar medium at rest, representing an explo- 
sion energy £sn ~ 10 51 erg and an ejected mass of 1.4 M Q . For 
a given distance the hydrogen density can then be derived from 
the known angular expansion velocity and size of the remnant, 
assumed to be given by the radio data of |Dickel et al. (1988)| and 
averaging these data over the azimuthal non-uniformities of the 
projected SNR shell. The use of such an average value for the an- 
gular velocity of the shock and the implied assumption of a uni- 
form circumstellar density is a necessary approximation within 
such a one-dimensional model which is meant to describe the 
overall physics of a point explosion. On the other hand, the sys- 
tematic errors which these assumptions introduce are difficult to 
estimate, especially in the transition between sweep-up and adi- 
abatic phase. BKV obtained the spectrum and the spatial distri- 
bution of CR in the remnant and the density of thermal gas. On 
this basis they then calculated the expected flux of non-thermal 
emission (Fig. [2]). To account for the uncertainties in the distance 
estimate this was done for a distance range from 3.4-7 kpc. The 
derived ambient density varies with the distance assumed and the 
numerical results show that for a distance d as low as 4.8 kpc the 
SNR has reached the Sedov phase. Therefore the predicted in- 
tegral hadronic gamma-ray flux roughly decreases with distance 
oc El^/d 1 , in agreement with the calculations shown in Fig. [2] 
Approximating the emission measure for free-free emission by 
EM ~ NnM sw , where M sw denotes the swept-up circumstellar 
mass, EM scales in the same way with £sn an d d as does the 
gamma-ray flux. 

To compare the given upper limit with the model predic- 
tion, the quantity F(> E) = E ■ F(> E) is determined. Here 
F(> E) is the upper limit on the integrated Flux above the en- 
ergy E. For E = 230 GeV the value for F is F{> 230 GeV) = 
3.4 X 10" 13 ergcirT 2 s~'. The resulting integrated upper limits 
are plotted in Fig. [2] for several energies in the range 0.23 < 
(E/l TeV) < 3.7. Note that for these values no upper bound for 
the integration is needed as the quantity F(> E) decreases with 
energy for spectral indices greater than F = 1 .0. 

Within the context of the model of BKV, the H.E.S.S. up- 
per limits rule out distances smaller than 6.4 kpc for £sn ^ 
10 5I erg and thus densities larger than 0.7 cm" 3 , and values of 
EM in excess of > 13M cirT 3 . The mean shock velocity is then 
~ 4000 km s _I , and the SNR is just in transition from the sweep- 
up phase to the Sedov phase. 

From SN explosion theory (see BKV, and references therein) 
a lower limit of £sn ~ 0.8 x 10 51 erg appears appropriate for 
type la SNe. Considering such a reduced explosion energy, the 
expected flux in gamma rays would be lower and therefore the 
H.E.S.S. upper limit would result in a reduced lower limit on the 
distance of d > 6.0 kpc. 

While in the BKV model the above value of EM that corre- 
sponds to the upper limit of H.E.S.S. agrees quite well with the 
overall number recently derived by Blair et al. (2007) from their 
measurements, the distances of 6.4 and 6.0 kpc differ signifi- 
cantly from the value adopted by these authors, whose distance 
estimate is within the errors smaller than 5.8 kpc (see Section[Tji. 
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Fig. 2. Comparison of the upper limits on F(E) with predictions 
of BKV. 



On the other hand, |Blair et al. (2007)| derived their distance value 
from an optical filament in the northwestern region which has 
the smallest expansion parameter found in the radio and X-ray 
observations all around the remnant. It is also interesting to note 
that the determination of the shock velocity from the Ha line 
broadening should involves a CR-modified shock, in contrast to 
the assumption of Sankrit et al. (2005) Efficient particle accel- 
eration in the SNR modifies the shock, whereby part of the gas 
compression - but only a very small part of the gas heating - 
occurs in a smooth precursor, in which the CR pressure gradi- 
ent slows down the incoming gas flow. This is followed by the 
so-called subshock (Dru ry & Volk 19 81) where most of the gas 
heating occurs (Berezhko & Ellis on 1999K The compression ra- 
tio cr s of the subshock is smaller than the overall shock com- 
pression ratio cr. Such a shock structure implies that the shock 
velocity corresponding to the downstream thermal motions is the 
subshock velocity V| ub = cr s /cr x V s , where V s is the total shock 
velocity. In other words, a higher overall shock velocity is re- 
quired to achieve the same gas heating if in addition CR are ac- 
celerated. Therefore the source distance derived from the width 
of the Ha line is cr^/cr < 1 times the true source distance if 
derived without particle acceleration. This may imply a substan- 
tial systematic error. In the BKV model for Kepler, in spherical 
symmetry it is cr s /cr « 0.4 for an assumed distance of 4.8 kpc, 
and still equal to 0.6 for d — 6.4 kpc. Therefore the nominal 
distance d — 4 kpc adopted by |Blair et al. (2007)| is equivalent 
to d — 6.6 kpc, if the northwestern region considered is indeed 
one where acceleration is efficient. If particle acceleration is not 
efficient in this region, then the correction factor is unity. Even 
a slight modification <r s /cr » 0.9 of the shock makes the source 
distances compatible. 

Independent of particle acceleration models one can use the 

H. E.S.S. upper limit also to constrain the content of energetic 
particles in the remnant. Using the limit on the flux in the range 
between 0.23 TeV and 12.8 TeV, an upper limit on the gamma- 
ray luminosity L y> max = 4-nFEd 2 can be estimated, where Fe is 
the integrated energy flux upper limit. In this range then L y max < 

I. x 10 32 • (d/kpc) 2 erg s _1 is derived. For power-law spectra, 
the (5-function approximation 



i \ cn EA E„ 
ME,)*- o- w -Up - 



(1) 



can be used to relate the spectra of pions (or gamma rays) to 
those of the primary protons ( jAharonian & Atoy an 2000); here 
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<p„ is the pion (or gamma-ray) production rate, n the gas den- 
sity, n p the number of protons and K„ is the mean fraction of 
the kinetic energy of the proton transferred to the secondary n°- 
meson per collision. The rest mass of the proton is neglected. 
For spectral indices F = 2 - 3, K n = 0.17 can be used to 
approximate the pion spectrum (Aharonia n & Atoyan 2000| , a 
similar value also applies for gamma-ray spectra. At high en- 
ergy, the proton-proton cross-section <x pp is only weakly en- 
ergy dependent and can be approximated with <x pp 40 mb 
dGaisser 19 90). The H.E.S.S. data probe proton energies in 
the range from about 1 to 100 TeV; using equation Q] and as- 
suming an index T = 2.0, a limit on the energy in pro- 
tons of £ p < 4.9 x 10 47 (c//kpc) 2 (7i/cirr 3 )~ erg can be de- 
rived. Extrapolating to the 1 GeV to 1 PeV range results in 
an upper limit on the total energy in accelerated protons of 

1.5 x 10 48 erg (d/kpc) 2 («/cirr 3 ) . Using d = 6.4 kpc and 
n = 0.7 crcT 3 (dashed-3-dotted line in Fig. [2j this results in £ p < 

8.6 x 10 49 erg, i.e. ~ 9% of the assumed energy £sn = er §- 
This is of the order of what is expected for an average cosmic- 
ray source in the form of a SNR. Assuming £sn = 10 51 erg and 
using the argument of BKV that d < 7 kpc, in agreement with 
the observational argument of |Reynoso & Goss (1999) the ex- 



pected gamma-ray flux should not be lower than the H.E.S.S. 
upper limit by more than a factor 2 as can be seen in Fig. [2] 

In another scenario the gamma-ray emission can be pro- 
duced via IC scattering by VHE electrons off ambient pho- 
tons mainly from the cosmic microwave background (CMB). 
The same electrons emit synchrotron X-ray radiation by be- 
ing deflected by magnetic fields in the SNR. The energy of the 
gamma-ray photons is coupled to that of the X-ray photons ac- 
cording to (Ex/lkeV) ~ 0.07 • (E 7 /l TeV) (B/10//G) in the 
case of the CMB as target photon field for the IC scattering. 
If the observed hard X-ray radiation (All en et al. 19991 1 , with 
a flux normalisation of 6.2 x 10~ 5 cm~ 2 s" 1 keV" 1 and a slope 
of -3.0 ± 0.2, is synchrotron radiation (with the correspond- 
ing energy flu?<Q fx) the energy flux in gamma rays is given 
by f r (E 7 )/fx(E x ) ~ £0.1(B/10^GT 2 ( I Aharonian et al. 1997D . 
The factor £ takes into account possible differences in the source 
sizes in X-ray and gamma-ray wavelengths. We assume here 

In principle one could try to use the above relations to ob- 
tain a lower limit on the magnetic field strength since the upper 
limit on the flux in particular constrains any IC component. For 
this purpose the energy flux from X-ray synchrotron emission 
at an energy corresponding to a given energy probed in VHE 
gamma rays has to be known either from measurements or from 
detailed modeling. The interval in Ex that corresponds to the 
observed gamma-ray energy interval 0.23 < E 7 /l TeV < 3.7 is 
~ (B/10juG) x (0.02 - 0.26) keV, whereas the energy interval in 
which the total non- thermal X-ray flux is known is 10 - 20keV 
d Allen et al. 199% . The X-ray instrument PC A on board RXTE, 
with which the underlying data were obtained, has no imaging 
capabilities and therefore the measured spectrum is the over- 
all spectrum of the field of view of the instrument (which is 
1°). Although it is expected that the measured photon flux is 
indeed from Kepler's SNR because of its position well above 
the plane, the X-ray flux has to be treated as an upper limit. 
Unfortunately there is no published analysis of the non-thermal 
flux from Chandra data covering the entire remnant. It is also not 
possible to unambiguously disentangle the non-thermal and the 



thermal contribution to the total spectrum measured by XMM- 
Newton (G. Cassam-Chenai', private communication). 

In the energy range around a few keV the extrapolation of 
the observed hard X-ray flux to lower energies involves consid- 
erable uncertainties. Nevertheless, in almost all scenarios the ex- 
trapolation of the power-law spectrum measured between 10 and 
20keV (with a spectral index of F = 3.0) should be an upper 
limit to the X-ray to UV flux. With this extrapolation an upper 
limit on the gamma-ray flux from IC scattering for a given mag- 
netic field can be calculated using the above formulas. 

For magnetic field values greater than 52 fiG the resulting 
predicted upper limit on the IC flux would be less than the mea- 
sured upper limit of Fe(3.7 TeV) < 2.91 x 10~ 13 erg cirT 2 s _1 . 

From Chandra measurements of thin X-ray filaments 
(Bambael al. 2005), whose thickness is interpreted as the syn- 
chrotron cooling length of the radiating electrons, the actual field 
strength is B ~ 300 fiG, following the arguments of BKV and 
|Parizot ef al. (2006) This field implies an IC gamma-ray energy 
flux of fy(Ey) < 1.4 X 10~ 15 erg cirT 2 s _1 which is two orders of 
magnitude below the measured upper limit. 

4. Conclusions 

Observations of Kepler's SNR with H.E.S.S. result in an up- 
per limit for the flux of VHE gamma rays from the SNR. In 
the context of an existing theoretical model (BKV) for the rem- 
nant, and assuming an ejected mass of 1 .4 M Q and an explosion 
energy of 10 erg in agreement with type la SN explosion mod- 
els, the lack of a detectable gamma ray flux implies a distance 
of at least 6.4 kpc, which is the same as the upper limit derived 
by|Reynoso & Goss (1999) from radio observations. Given that 



the gamma-ray flux effectively scales with £^ N , a significantly 
higher explosion energy is excluded; a theoretically acceptable 
lower explosion energy of 0.8 x 10 erg would lower the distance 
limit to 6 kpc. 

Assuming a purely hadronic scenario, a standard type la SN 
explosion, and using 6.4 kpc as a lower limit for the distance, the 
H.E.S.S. upper limit implies that the total energy in accelerated 
protons is less than 8.6 x 10 49 erg. 

In a synchrotron/IC scenario no strong constraints on the 
magnetic field can be obtained. 
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